Phonon and plasmon excitation in inelastic electron tunneling spectroscopy of 
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The inelastic electron tunneling spectrum (IETS) of highly oriented pyrolitic graphite (HOPG) has 
been measured with scanning tunneling spectroscopy (STS) at 6K. The observed spectral features 
are in very good agreement with the vibrational density of states (vDOS) of graphite calculated 
from first principles. We discuss the enhancement of certain phonon modes by phonon-assisted 
tunneling in STS based on the restrictions imposed by the electronic structure of graphite. We also 
demonstrate for the first time the local excitation of surface-plasmons in IETS which are detected 
at an energy of 40 meV. 

PACS numbers: 68.35.Ja, 68.37.Ef, 63.20.-e 



Initially applied to measure inelastic tunneling cur- 
rents across metal- insulator-adsorbate-metal junctions 1 , 
Inelastic Electron Tunneling Spectroscopy (IETS) has in 
the recent past been revived and its concept extended to 
Scanning Tunneling Microscopy (STM-IETS)^. Super- 
imposed on the elastic tunneling current, the inelastic 
tunneling process opens up alternative tunneling chan- 
nels increasing the total conductance at the onset of its 
excitation. Very recently, elegant experiments (e.g i 2 i 3 i 4 i 5 ) 
have been performed in detecting local inelastic electron 
energy-loss spectra of individual molecule adsorbed on 
metallic substrates. The observed features in the tun- 
neling spectra are found to correspond to the vibra- 
tional modes of the adsorbed species. Vibrations of single 
molecules as large as Ceo have have been detected^. De- 
spite the increasing number of publications of STM-IETS 
experiments on single molecules adsorbed on metals, to 
our knowledge, only two studies were reported on collec- 
tive vibrational excitations of surfaces^. In both these 
studies graphite was used as a substrate and several spec- 
tral features remained unexplained or were attributed to 
vibrational excitation in the tungsten tip. However, these 
tip effects have not been detected in many other STM- 
IETS experiments. In order to resolve the pure inelastic 
response of a graphite surface, we have performed STM- 
IETS experiments on a highly oriented pyrolitic graphite 
crystal (HOPG) at cryogenic temperature. In our experi- 
ments the whole vibrational spectrum of HOPG from the 
rigid layer shear mode at « 5meV to the optical modes at 
« 200 meV is clearly observed without being perturbed 
by inelastic effects from the tip. The phonon modes of 
graphite can be clearly identified by the direct compar- 
ison of experimental data with first principles density 
function theory (DFT) calculations. We report for the 
first time that also surface-plasmon losses can also be 
observed directly with STM-IETS. 

Freshly cleaved HOPG has been introduced into the 



vacuum system and shortly annealed at 500K. The ex- 
periments were performed at 6K with a home built UHV- 
STM using a tungsten tip. The inelastic tunneling spec- 
trum of HOPG was measured over a topographically 
clean and flat terrace with a lock-in technique. In this 
technique a sinusoidal signal with amplitude AV is su- 
perimposed to the applied sample bias. In IETS the sec- 
ond harmonic component of the modulated signal, which 
is proportional to the second derivative vs. bias of the 
tunneling current, d 2 I/dV 2 , is recorded. In this way the 
inelastic loss features in the conductance are transformed 
into peaks at the onset of the loss modes for positive 
sample bias and into minima for negative sample bias, 
respectively. A typical example for HOPG is shown in 
fig.l. 




FIG. 1: Inelastic electron tunneling spectrum of HOPG mea- 
sured by STM. The phonon features (arrows) in the d 2 I/dV 2 
are symmetrically located with respect to the zero bias posi- 
tion. 



We obseve phonon features that are symmetrically lo- 
cated with respect to zero bias as expected for the inelas- 
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tic tunneling process. Possibly due to a particular tip-end 
configuration teh phonon on the negative bias side are 
more pronounced. In Fig. 2a we present a high-resolution 
spectrum of the d 2 I/dV 2 signal of HOPG taken at nega- 
tive bias. The data was obtained with a bias modulation 
of AV = 10 meV RMS at a frequency of 2.7kHz and 
inverted to render the phonon excitations as peaks to 
facilitate the comparison to theory. Since the electron 
density of states of HOPG is flat in the energy range 
close to the Fermi level, the features in the spectrum in 
Fig. 1 are exclusively due to inelastic processes. 
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FIG. 2: a) Inelastic electron tunneling spectrum of HOPG 
measured by STM. Inelastic excitations show up as peaks 
in the d 2 I/dV 2 signal. b) vibrational density of states 
calculated by DFT-LDA. c) Phonon dispersion relation of 
graphite. The phonon branches are specified as following: 
out of plane acoustic (ZA), acoustic shear (SH), longitudinal 
acoustic (LA), out of plane optical (ZO), optical shear (SH*), 
and longitudinal optical (LO). The inset shows the first Bril- 
louin zone of graphite with the high symmetry points. A 
detailed description of the observed phonon modes is given in 
the text. Furthermore, a plasmon excitations is identified at 
40 meV. 

Along with the measured inelastic spectrum, we 
present in Fig. 2b and c the calculated vibrational den- 
sity of states (vDOS) and the phonon dispersion relation 
of graphite along the A-F-M-K directions. Graphite is a 
layered structure with strong carbon bonds in the plane 
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TABLE I: Observed features in the STM-IET spectrum of 
HOPG in comparison with vibrational modes observed in ex- 
periment (Neutron scattering^, HREEL S 12 ' 13 and double- 
resonant Raman scattering—) and with first principles DFT- 
LDA calculations. 



and weak Van der Waals-like forces between the planes. 
As a result of this anisotropy, the phonon spectrum of 
graphite covers a very wide range. We have performed 
first principles calculations using density-functional the- 
ory (DFT) in the local density approximation (LDA) 9 . 
Aim of this calculation is to resolve some of the devia- 
tions in the vDOS calculated by the different recent force- 
constant parametrizations for graphiteiSiiiiiii^ and to 
have a solid foundation for the assignment of the mea- 
sured peaks. The first principles calculation has been 
shoAvi^M^I— to yield very close agreement with phonon 
dispersion relations measured by high resolution elec- 
tron energy loss spectroscopy (HREELS) along the T-M 
direction^ and along the T-K direction^. Only for the 
acoustic shear-mode (SH) it was observed^ that at the 
M-point, the HREELS data of Refill yields an energy 
which is higher by 25 meV than the ab-initio value. A 
similar underestimation of of the shear-mode frequency 
might occur as well at the K-point where a mode with 
energy 134 meV was observed by double-resonance Ra- 
man scattering^. Our calculated vDOS is also at vari- 
ance with the HREELS data at Mi£ but is in very good 
agreement with the vDOS measured by neutron scatter- 
ing on a powdered sample of graphite 2 » where a separate 
SH(M)-feature was not observed. However, since the SH 
K-point mode contributes only weakly to the calculated 
DOS the comparison to the experimental vDOS in Ref. 
l2ll cannot answer the question whether the mode comes 
out too soft in the calculations. 

The observed vibrational spectral features, their com- 
parison with values obtained with other experimental 
techniques and their assignment to phonon modes at 
the various points of the surface Brillouin zone (SBZ) 
are summarized in Table 1. Most features coincide with 
peaks in the calculated vDOS which arise from the flat 
dispersion relation around the high-symmetry points. In 
particular, the higher peaks due to phonon modes around 
the M-point are clearly mapped. The first peak at 
16 meV is in good agreement with the calculated energy 
of the out of plane acoustic (ZA) mode at T where two 
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neighboring planes are oscillating out-of-phase. While in 
the vDOS and in neutron scattering^ this peak is only 
weak, in the IETS experiment it is stronlgy pronounced 
which can be explained by the direction of oscillation 
which is orthogonal to the surface plane. The weak shoul- 
der at 6 meV corresponds to the low frequency horizontal 
mode at T where two neighboring planes are oscillating 
out-of-phase parallel to each other. The prominent fea- 
ture at 40 meV is not due to a vibrational feature but is 
assigned to a plasmon excitation (see below). In RefQ, 
a phonon feature at 39 meV was observed which was as- 
sociated with a double-phonon excitation in the surface 
(16 meV) and in the tungsten tip (23 meV). However, in 
our experiment the tip-phonon peak is clearly not present 
and all the observed features can therefore be assigned 
to surface excitation. The vibrational feature observed at 
58 meV can be assigned to the ZA at the M-point. The 
optical out-of-plane (ZO) mode and the shear acoustic 
mode SH (also called transverse acoustic mode) display 
a saddle in the dispersion relation around the M-point 
adding up to a very high peak in the vDOS at 78 meV. 
This explains the peak at 81 meV in the measurement. 
However, as noted above, HREELS datai^ suggests the 
energy of the SH mode at M to be about 100 meV. This 
gives rise to the alternative explanation that the peak at 
81 meV is only due to the ZO mode while the shoulder 
at 100 meV may be a signal of the SH mode at the M- 
point. The ZO mode at T which gives rise to a peak at 
108 meV in the theoretical vDOS seems to be missing in 
the spectrum. Instead, the spectrum reveals actually a 
dip (dl 2 /dV 2 < 0) at this energy which will be discussed 
below. At 124 and 150 meV two low-intensity phonon 
modes, the SH, and the longitudinal acoustic (LA) modes 
at the K-point are predicted, respectively. Experimen- 
tally we have observed a large asymmetric inelastic fea- 
ture at about 137 meV. This disagreement between ex- 
periment and theory might point to an underestimation 
of the SH mode also at the K point. If, as could be the 
case for the SH mode at Mi£, the real energy is higher by 
10 to 20 meV, the peak at 137 meV can be assigned to 
the SH mode at K which corresponds to the above men- 
tioned observation of a K-point phonon at 134 meV— 
by double resonance Raman spectroscopy. In the highest 
energy range of the spectrum two loss-features protrude 
the background signal at 180 and 200 meV, respectively. 
The first is related to the optical shear mode (SH*, also 
called transverse optical mode) around the M-point. The 
neighboring feature at 200 meV corresponds to the lon- 
gitudinal optical (LO) branch which displays a strong 
overbending in the direction T-M and T-K giving rise to 
a peak in the vDOS. 

In general, the measured IET spectrum follows very 
closely the vibrational density of states with two no- 
table exceptions which will be discussed in the follow- 
ing. These exceptions are related to the question if all 
phonons can be excited or if the geometry of the tunnel- 
ing system leads to a restriction to phonons in a certain 
region of the Brillouin zone, as it is the case in phonon- 



assisted electronic tunneling through narrow p-n junc- 
tions of semiconductors. In e.g. Ge, the valence-band 
maximum is centered at T while the conduction-band 
minimum is located at the L point. This allows only vi- 
brational modes at L to be excited by a phonon-assisted 
tunneling process^. The electronic structure of graphite 
is also quite unique being a semimetal with a band cross- 
ing only at the K-point of the SBZ. Therefore the elastic 
current is carried by electrons that tunnel with a rela- 
tive large fcy = K (1.7 A -1 ). In an STM geometry the 
expansion of the initial state, i.e., the tip wave function 
in terms of surface plane-waves^, leads to a parallel mo- 
memtum distribution centered at feii = 0. Consequently 
the electrons of the elastic current have a low transmis- 
sion probability. Electrons with any other k» can only 
tunnel by exchange of finite crystal momentum with a 
phonon. The phonon-assisted tunneling is relevant to all 
the dominant features observed besides the T phonons 
at 16 and 200 meV. An enhanced tunneling contribution 
should especially occur for electrons with initial paral- 
lel momentum fcy = that tunnel to the Fermi level at 
the K-point by exciting (or absorbing) a K-point phonon. 
This explains why the peak at 137 meV which is probably 
due to a K-point phonon excitation shows up so domi- 
nantly in the spectrum even though the vDOS of states 
only displays weak features for K-point phonons. 

The second deviation of the measured spectrum from 
the vDOS of graphite is the dip at 111 meV where the 
vDOS displays a peak due to the ZO phonon at T. The 
theory for inelastic tunneling spectroscop}i22i2£ distin- 
guishes between an inelastic contribution (real phonon 
emission or absorption) and an elastic contribution (vir- 
tual phonon emission with subsequent reabsorption or 
phonon absorption with subsequent re-emission). In 
some cases, the elastic contribution leads to a decrease 
of the tunneling conductance (dl/dV) due to back- 
scattering of the electron^. This effect has been mea- 
sured recently as dips in the dI 2 /dV 2 curve of O2 on a 
silver surface^. We believe that the dip at 111 meV is 
due to a similar elastic contribution at the onset of vDOS 
singularity corresponding to the ZO(r) phonon branch. 
In the light of the above discussion of phonon-assisted 
tunneling on graphite, it is also not surprising that this 
effect is observed for phonons at T. The only matching 
pair for a virtual phonon scattering process are electrons 
initially at ku = K and phonons at T. All other combina- 
tions are forbidden due to the unavailability of intermedi- 
ate electronic states. Furthermore, the steady reduction 
of -d 2 I/dV 2 from the ZO (M) peak at 81 meV, in con- 
trast to the vDOS, shows that symmetry of the phonon 
displacement (the ZO branch) should play an important 
role in the supression of the IET curreirfcSI. We note that 
one of the degenerate modes of the ZO branch at V is the 
only r phonon mode that has A2 U symmetry. 

In addition to the predicted phonon modes, the mea- 
sured spectrum contains an energy loss feature at 40 
meV, which cannot be explained by phonon excitations. 
Other than the previous features, this loss feature is 
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of purely electronic origin and corresponds to an out- 
of-plane oscillation of the electrons. A tunneling elec- 
tron traversing a vacuum-metal interface has a signif- 
icant probability of interacting with collective electron 
oscillations. The overlapping of the tt and it* band at 
the K point of the SBZ yields a low density of free car- 
riers (i.e. electrons and holes) near the Fermi level and 
allows low-energy plasmon excitation. A strong tempera- 
ture dependence of the excitation energy of this plasmon 
mode was observed in HREELS 28 . The excitation en- 
ergy shifts from 63 meV to 45 meV when going from 
room temperature to 150 K, respectively. This effect has 
been attributed to the low density of states at the Fermi 
level in conjunction with the thermal excitations of car- 
riers. Thermal excitations from the valence into the con- 
duction band can radically alter the electron population 
at the Fermi level and yields the strong dependence of 
surface plasmon energy and line width on temperature. 
The temperature dependence of this excitation energy 
was also calculated in Ref. |28[ Within the random phase 
approximation the plasmon energy excitation was pre- 
dicted to be lowered to 36 meV at OK in rather good 
agreement with our experimental value. The low-energy 
plasmon was also observed in transmission electron en- 
ergy loss spectroscopy, which also detected a second plas- 
mon loss, assigned to an excitation polarized parallel to 
the layer at 128 meV— . The shoulder observed around 



125 meV might be assigned to this plasmon excitation. 

In conclusion, the IET spectrum of HOPG has been 
measured by low temperature scanning tunneling spec- 
troscopy and is compared to an ab-initio calculation of 
the vibrational density of states in graphite. The ob- 
served vibrational loss features can be identified with 
peaks in the vDOS which are due to the flat dispersion of 
the phonon branches around the high-symmetry points 
of the Brillouin zone. Additionally, a prominent plasmon 
feature at 40meV has been identified in the inelastic tun- 
neling current measured with STM. In the past, the local 
capability of STM-IETS has been proven to be a valuable 
tool to access the vibrational properties of adsorbates at 
surfaces. In this work we have shown that STM-IETS re- 
solves the vibrational and collective electronic excitations 
of graphite with high accuracy. This is an important step 
for the application of the method to confined geometries 
and individual carbon nanostructures. 
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